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Smart Pilot Assignment for Massive MIMO 

Xudong Zhu, Zhaocheng Wang, Linglong Dai, and Chen Qian 


Abstract —A massive multiple-input multiple-output (MIMO) 
system, which utilizes a large number of antennas at the base 
station (BS) to serve multiple users, suffers from pilot contamina¬ 
tion due to inter-cell interference. A smart pilot assignment (SPA) 
scheme is proposed in this letter to improve the performance of 
users with severe pilot contamination. Specifically, by exploiting 
the large-scale characteristics of fading channels, the BS firstly 
measures the inter-cell Interference of each pilot sequence caused 
by the users with the same pilot sequence in other adjacent 
cells. Then, in contrast to the conventional schemes which assign 
the pilot sequences to the users randomly, the proposed SPA 
method assigns the pilot sequence with the smallest inter-cell 
interference to the user having the worst channel quality in a 
sequential way to improve its performance. Simulation results 
verify the performance gain of the proposed scheme in typical 
massive MIMO systems. 

Index Terms —Massive MIMO, pilot contamination, inter-cell 
interference, pilot assignment, channel quality. 


I. Introduction 

Massive multiple-input multiple-output (MIMO) has been 
recently investigated to meet the exponential increase of mo¬ 
bile traffic in wireless systems m, El, whereby a base station 
(BS) equipped with a large number of antennas serves multiple 
users simultaneously. Asymptotic analysis based on random 
matrix theory ^ demonstrates that the intra-cell interference 
and the uncorrelated noise can be effectively eliminated when 
the number of BS antennas goes to infinity. However, pilot 
contamination caused by the inter-cell interference originated 
from the reuse of the same pilot group in adjacent cells, does 
not vanish as the increase of BS antennas, and it becomes the 
performance bottleneck of massive MIMO systems fTl, ||2|. 

The issue of pilot contamination has been widely studied 
in the literature 0-0- The time-shifted pilot scheme is an 
effective solution by using asynchronous transmission among 
adjacent cells 0, but leads to the mutual interferences be¬ 
tween data and pilot. A greedy pilot assignment algorithm 
0 can mitigate the pilot contamination by exploiting the 
statistical channel covariance information, but suffers from 
high computational complexity of iteratively minimizing the 
mean degree of spatial orthogonality. Pilot contamination 
precoding 0 can mitigate the inter-cell interference by multi¬ 
cell joint processing, but suffers from spectral efficiency loss 
due to high overhead required by information exchange. In 
addition, a blind method based on subspace partitioning 0 
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Fig. 1. A typical multi-cell multi-user massive MIMO systems. 

is able to reduce the inter-cell interference when the channel 
vectors from different users are orthogonal. All those solutions 
assign the available pilot sequences to different users randomly 
without considering their different channel qualities, and also 
ignore the fact that the severity of pilot contamination varies 
among different pilot sequences. 

In this letter, a smart pilot assignment (SPA) scheme is 
proposed to enhance the performance of users with severe 
pilot contamination. We consider the pilot assignment problem 
for a target cell, which is surrounded by other adjacent cells. 
Unlike the conventional schemes which assign the available 
pilot sequences to the users in a random way, our proposed 
SPA scheme aims to maximize the minimum uplink signal-to- 
interference-plus-noise-ratio (SINR) of all users in the target 
cell. Specifically, by exploiting the large-scale characteristics 
of fading channels, the BS firstly measures the inter-cell 
interference of each pilot sequence caused by the users with 
the same pilot sequence in other adjacent cells. After that, the 
channel qualities from different users in the target cell to the 
BS can be detected, which usually differ from one user to 
another. The proposed SPA method assigns the pilot sequence 
with the smallest inter-cell interference to the user having 
the worst channel quality in a sequential way until all users 
have been assigned by their corresponding pilot sequences. 
Simulation results verify the effectiveness and performance 
gain of the proposed SPA scheme in typical massive MIMO 
systems. 



II. System Model 

As shown in Fig. [T] we consider a multi-cell multi-user 
massive MIMO system composed of L hexagonal cells, and 
each cell is consisted of a BS with M antennas and K (K -C 
M) single-antenna users Q, 0. The channel vector £ 
qMxi ffojjj fjjg ^sgj. jjj fjjg j.fjj ggjj (o the BS in the i-th 
cell can be modeled as 
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hijfc — Sij k's/ l^ij k 1 (1) 

where Pijk denotes the large-scale fading coefficients which 
change slowly and can be easily tracked 0-0, and gijk ~ 
CAf{0, \m) denotes the small-scale fading vectors. By consid¬ 
ering the typical protocol of time-division duplexing (TDD) 
in massive MIMO systems, we adopt the widely used time 
block fading model, whereby the channel vector remains 
constant during the coherence interval 0, 0. 
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Assuming the pilot sequences $ = 02: ’ ’ ’ ^ 

with length of r used in one cell are orthogonal, i.e., 
= Ik, and the same pilot group is reused in other cells 
due to the limited pilot resource m. The conventional pilot 
assignment methods usually assign the pilot sequence 0^, to 
the fc-th user without considering the different channel qual¬ 
ities among users m, Thus, the received pilot sequence 

yP g cMxr jjjg gs 

in the z-th cell can be represented as 

L K 

+N?, (2) 

2 = 1 fc=l 

where pp denotes the pilot transmission power, and S 
qMxt (jenotes the additive Gaussian white noise (AWGN) ma¬ 
trix with entries being independent and identically distributed 
(i.i.d.) Gaussian random variables with zero-mean and variance 
ajf. Similarly, the received user data y“ G at the BS in 

the z-th cell can be represented as 


L K 



2 = 1 fc=l 


where Xj/. denotes the symbol from the fe-th user in the j-th 
cell with E{|a;“^,p} = 1, pu denotes the uplink data trans¬ 
mission power, and n“ G denotes the AWGN vector 

with £’{n“(ni)^} = By correlating the received pilot 

sequence with pilot sequence 0^,, the channel estimate of 
the fc-th user in the z-th cell can be represented as 

1 ^ 

hufc = —YP0f = ^ + v,fc, (4) 

VPp j=i 

where = --^N^0^ denotes the equivalent noise. It is 
clear that the channel estimate of the fc-th user in the z-th 
cell, hiik, is a linear combination of the channels hy^, j = 
1,2, • • • ,L, of the users with the same pilot sequence in all 
cells, which is referred to as pilot contamination m, m. 

By adopting the matched-filter (MF) detector based on the 
channel estimate result huk, the detected symbol for the /c-th 
user in the z-th cell can be represented as 

^ik ^iiky i 

= (Eh.2fc +Wfe)^(v^E E +<) 

2 = 1 2 = lfe'=l 

= + E ^fjkKkX%) + e“fc> (5) 

2/* 


where denotes the intra-cell interference and uncorrelated 
noise, which can be significantly reduced by increasing the 
number of BS antennas Q. Then, the uplink SINK of the 
/c-th user in the z-th cell can be calculated as 


SINR“fc = 


Ihffeh, 
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( 6 ) 

Thus, the corresponding average uplink capacity of this user 
can be calculated as C“;. = E{log 2 (l -f SINR“j,)}. It is clear 
that thermal noise and small-scale fading effects could be 
averaged out as M grows to infinity. However, the average 


uplink capacity is limited by the pilot contamination and 
cannot be improved by the increase of either p^ or pp. 

III. Proposed Scheme 

In this section, the pilot assignment for a target cell is firstly 
formulated as an optimization problem. Then, the SPA scheme 
is proposed to approach the optimization solution in a greedy 
way. The performance analysis is also provided to verify the 
effectiveness of the proposed SPA scheme. 


A. Problem Formulation 


We consider the pilot assignment for a specific cell, i.e., 
the z-th cell as the target cell, and the pilot assignment for 
other cells are independently managed by their corresponding 
BSs. For this target cell, the number of different kinds of pilot 
assignments between K users [Ui,U 2 ,--- ,Uk] and K pilot 
sequences [0^^, 02, • • ■ , 0/f] is huge, i.e., P{K, K) = Kl. In 
the conventional pilot assignment schemes Q, El, the pilot 
sequence 0^. is assigned to the fc-th user Uk randomly. Based 
on the possible cooperation among cells, such as coordinated 
multiple points (CoMP) in LTE-A systems, a specific BS is 
able to acquire the pilot assignment knowledge of other cells. 
Without loss of generality, for other cells around the target 
cell, we assume that the random pilot assignment is used. 

In massive MIMO systems, since the poor users with severe 
pilot contamination are the performance bottleneck HI, 121. we 
aim to maximize the minimum uplink SINR of all K users 
in the target cell, which can be formulated as the following 
optimization problem 


V 


: max min 
{.F.} Vfe 



(7) 


where : s = 1, • • • ,7^!} denotes all possible K\ kinds 
of pilot assignments, e.g., Ps = Ifs’fsG'' , f^] denotes the 
s-th assignment, and assuming that the pilot sequence 0^, is 
assigned to the fc-th user Uk in all other cells. 

However, it seems impossible to solve this optimization 
problem V due to the fact that we cannot obtain accurate 
channel estimate under pilot contamination as shown in (01. 
Fortunately, the limit of the uplink SINR can be represented 
by the large-scale fading coefficients Pijk as shown in (|6]l. 
As shown in a-in, the large-scale fading coefficients Pijk 
change slowly and can be easily tracked by the BSs. Thus, 
the optimization problem V can be approached by 


V 


M—¥00 




: max min 
Vfc 



Ej/i P 


2 

ijk 


( 8 ) 


The proposed SPA scheme aims to solve this optimization 
problem V in a greedy way, which will be addressed in detail 
in the next subsection. 


B. Pilot Assignment 

The most direct way to solve the optimization problem V' 
is the exhaustive search, which tries all possible assignments 
and chooses the best one. However, the number of all pilot 
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assignments is as huge as K\, which leads to high computa¬ 
tional complexity. In this letter, the optimization problem V' 
is solved in a greedy way with low complexity. 

For the target cell, we define a series of parameters 
to quantify the channel quality of K users as 

ak=lilk, k = l, 2 ,--- ,K. ( 9 ) 

For the K pilot sequences 02; ‘ ‘ define an¬ 

other series of parameters {'yk]k=i quantify the inter-cell 
interference of each pilot sequence caused by the users with 
same pilot sequence in other adjacent cells as 

fc=l,2,... ,/F, (10) 

which varies among K pilot sequences. 

For a specific assignment Ts = [fl, ‘' i /i^]> the pilot 
sequence 0^, is assigned to the user Ufk, whereby the limit 
of uplink SINR of the user Ufk, i.e., SINR“j, —?► afkjjk, is 
decided by two aspects; 1) the channel quality ap of the 
user Ufk \ 2) the inter-cell interference 7^ caused by the users 
with the same pilot sequence 0 ^ in adjacent cells. In order to 
maximize the minimum uplink SINR of all users in the target 
cell, we need to avoid the pilot sequence with great inter-cell 
interference assigned to the user having bad channel quality, 
which leads to relatively low uplink SINR. 

Based on this motivation, the SPA scheme is proposed, 
which assigns the pilot sequence with the smallest inter-cell 
interference to the user having the worst channel quality in a 
sequential way. Mathematically, as shown in Fig. | 2 ] we first 
sort the K pilot sequences according to the severity of inter¬ 
cell interference in descending order, i.e., 

Tp : [4>fi , 0 ,---, 0 jk ], ( 11 ) 

whereby this permutation Tp of the K pilot sequences satisfies 
7/1 > 7/2 > • • • > 7/«; > 0 . Similarly, we sort the K users 
according to their channel qualities in descending order, i.e., 

( 12 ) 

whereby this permutation J^g of the K users satisfies ct/i > 
a/2 > ••• > a/K > 0 . Then the pilot sequence (pfk 
will be assigned to the corresponding user Ufk. Apparently, 
the computational complexity of the proposed SPA scheme 
comes from the sorting process, which is only 0 {K log K) 
and negligible compared with C 9 (Ar!) required by exhaustive 
search. 

In order to apply the proposed greedy solution to the whole 
system with L cells, a sequential iterative scheme can be 
utilized. Specifically, all L cells solve their own optimization 
problems in a sequential way, and then this sequential proce¬ 
dure is carried out iteratively until the convergence is achieved. 

C. Performance Analysis 

In this subsection, we prove that the pilot assignment 
greedily generated by the proposed SPA scheme is one of the 
solutions to the optimization problem V. 

Without loss of generality, we fix the permutation of K 
pilot sequences as Fp addressed in (fTTT i. and assume that 
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Fig. 2. Intuitive description and proof of the proposed SPA scheme. 


there is another permutation J7/ of K users, which has larger 
minimum uplink SINR of K users than that of J7, i.e.. 


OLjk 

^q' = min —— > /i„ = min —( 13 ) 
Vfc 7jfe Vfe yjk 

Now we prove that such permutation Tg' does not exist 
and the intuitive proof is shown in Fig. | 2 ] More specifically, 
we firstly assume that the user U .kg has the minimum uplink 

Jq 

SINR in the permutation i.e., = a,kg/y .kg. Then, 

Jq Jp 

according to the permutation J^qf, the pilot sequence cp.ko will 

Jp 

be assigned to the user U.kg . Comparing a .kg with a .kg , there 
Jq' N' ■'« 

will be two results: 


1 ) a.kg < a.kg: Due to the fact that the permutation of 

Jq' J'J 

the pilot sequences J^p is fixed, we have 


Pq' < 


Oi. rkQ 
^ a' 


OL ^A:q 


< 


— 


( 14 ) 


which contradicts with the assumption in (fT 3 ]) . 

2 ) a .kg > a .kg : Due to the sorting process in (fTSl i. it is 
Jq' U 

clear that 


fg^ Jq°-^}- ( 15 ) 

Then, we consider the pilot assignment of the pilot 
sequences {cf>fi, 0x2, • • • , cj).kg-i}. Since one user cho- 
sen from {Ufi,Uf2, ■ ■ ■ ^U^kg-i} has been assigned 
with the pilot sequence 0,^0, there must be an user 

Jp 

U.kq (ki > ko) assigned with a pilot sequence (f.kg 

Jq Jp 

(k2 < ko). As a result, we have 

Qi CX pkQ 

Mg' < ^ ( 16 ) 

Jp Jp Jp 

which contradicts with the assumption in ( [ 13 ] ). 

It is concluded that the permutation IFg generated by the 
proposed SPA scheme is one of the solutions to the opti¬ 
mization problem V'. Note that there may be more than one 
optimal solution, for example, there are two users with the 
same channel quality, i.e., Ukq = akgjki k2. However, the 
exchange of the pilot sequences assigned to these two users 
makes no difference on the minimum uplink SINR. 

It should be pointed out that the solution to the optimization 
problem V is also the solution to the optimization problem V 
when the number of antennas at the BS grows to infinity, i.e., 
M ^ 00. With a finite number of BS antennas in practical 
systems, the solution to V' can approach the solution to V, 
which will be verified by simulation results. 
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Fig. 3. Simulation results: (a) The CDF of the minimum uplink SINR among all K users in the target cell; (b) The average uplink capacity of the user with 
the minimum uplink SINR among in the target cell; (c) The convergence of the average minimum uplink SINR in each cell. 


IV. Numerical Results 

In this section, we investigate the performance of the 
proposed SPA scheme through Monte-Carlo simulations. A 
typical hexagonal cellular network with L cells is considered, 
where each cell has K users with single-antenna and a BS 
with M antennas m, m. The center cell surrounded by other 
cells is considered as the target cell. The system parameters 
are summarized in Table I. As addressed in m, the large-scale 
fading coefficient Pijk can be modeled as 

Pijk = Zijkl{rijklRT, (17) 

where Zijk represents the shadow fading and possesses a log¬ 
normal distribution (i.e., 10 log]^Q(zi_,fe) is Gaussian distributed 
with zero mean and the standard deviation of (Tshadow), fijk is 
the distance between the fc-th user in the j-th cell and the BS 
in the i-th cell, and R is the cell radius. 


TABLE I 

Simulation Parameters 


Number of cells L 

7 

Number of BS antennas M 

8 < M < 512 

Number of users in each cell K 

K = 8 

Cell radius R 

500 m 

Cell edge SNR 

20 dB 

Average transmit power at users pp,pu 

0 dBm 

Path loss exponent a. 

3 

Log normal shadowing fading (Jshadow 

8 dB 


Fig. [3 (a) plots the cumulative distribution function (CDF) 
curve of the minimum uplink SINR among all K users in the 
target cell. The conventional schemes assign the pilot sequence 
to the user Uk (HI. 121. the optimal solution is the 
solution to the optimization problem V obtained by exhaustive 
search. When M = 32 is considered, it is evident that the 
proposed SPA scheme outperforms the conventional schemes 
by about 2 dB, and the optimal solution is the best with another 
gain about 0.5 dB. When the number of BS antennas increases, 
i.e., M = 512 is considered, we can find that the performance 
of the proposed SPA scheme is almost the same as that of the 
optimal solution. 

Fig. |3](b) shows the average uplink capacity of the user with 
the minimum uplink SINR among all K users in the target 
cell. When the typical parameter M = 128 is considered, the 
average uplink capacity of the user with the minimum uplink 


SINR in the proposed SPA scheme is larger than that of the 
conventional schemes by about 0.6 bps/Hz. Moreover, when 
M grows large, the performance of the proposed SPA scheme 
approaches that of the optimal solution. 

Fig. [3 (c) shows the convergence of the average minimum 
uplink SINR in each cell when the sequential iterative scheme 
is applied to the whole system. When M = 32 is considered, 
it takes about 5 iterations to converge for all L cells, and the 
required number of iterations becomes smaller if the number 
of BS antennas increases, e.g., M = 128 or M = 512 as 
illustrated in Fig. [3](c). 

V. Conclusions 

In this letter, we have proposed a smart pilot assignment 
scheme to improve the minimum uplink SINR of all users 
within the target cell in massive MIMO systems. By exploiting 
the large-scale characteristics of fading channels, the proposed 
SPA scheme assigns the pilot sequence with the smallest 
inter-cell interference to the user having the worst channel 
quality in a sequential way. Theoretical analysis proves that 
the pilot assignment generated by the SPA scheme is the 
solution to the simplified optimization problem V', which can 
approach the original optimization problem V when M grows 
to infinity. Simulation results demonstrate that for the typical 
configuration of 128 BS antennas in multi-user massive MIMO 
systems, the proposed SPA scheme is able to improve the 
minimum uplink SINR of users by about 2 dB. 
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